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EXECUTIVE  SUMMARY 


Induction  motors  are  used  worldwide  in  many  residential,  commercial,  industrial 
and  utility  applications.  Recently,  they  are  the  preferred  choice  among  the  industrial 
motors  due  to  the  modern  power  electronics  that  improve  their  speed  control.  Vitally 
important  for  the  speed  control  of  a  motor  is  the  accurate  estimation  of  the  magnetic  flux 
and  the  electromagnetic  torque.  Knowing  the  electromagnetic  torque  of  a  motor,  we  are 
able  to  control  it  and  thus  monitor  the  speed  faster  and  more  stably. 

The  purpose  of  this  study  is  to  design  an  electromagnetic  torque  estimator  of  a 
three-phase  induction  motor,  without  the  use  of  any  sensor.  This  goal  is  being  achieved 
by  building  the  mathematical  model  of  the  induction  motor  as  well  as  deriving  the 
electromagnetic  torque  equation  by  means  of  the  proper  transformations  and 
mathematical  tools.  The  analysis  of  the  designed  estimator  is  done  through 
Simulink/Matlab,  where  the  transient  and  steady  state  torque  of  the  motor  is  detennined. 
Furthermore,  the  experimental  test  of  the  estimator,  using  a  real  motor,  power  electronics 
and  the  Field  Programmable  Gate  Array,  proves  its  efficiency. 

The  electromagnetic  torque  equation  of  the  motor  in  the  arbitrary  (qdO)  reference 
frame  becomes: 

Te=~[iqsKs-ids\s)  (a) 

In  the  equation  above: 

•  P  stands  for  the  number  of  the  motor  poles. 

•  The  quantities  i  and  ids  stand  for  the  q-axis  and  d-axis  currents 
respectfully. 

•  A  andT/(  stand  for  the  magnetic  fluxes  on  the  q-axis  and  d-axis. 

Equation  (a)  can  be  transformed  into  the  stationary  reference  frame  where  the 

speed  ( a  =  d0  )  is  zero  as: 
dt 
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T.  =  || (rA*  -iM* )  =||{r  (| (v*  (J (v«-  -I-/1)*)}  (b) 


The  torque  equation  (b)  depends  on  the  transfonned  stator  voltages  (v^andv*), 
the  stator  currents  ( iqs  and  ids )  and  the  resistance  of  stator  phases.  Accordingly,  the 
formulation  of  the  torque  equation  is  based  on  equation  (b),  and  in  the  Simulink/Matlab 
environment,  it  can  be  designed  as  shown  in  Figure  1 : 
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Figure  1.  Torque  Estimator  in  Simulink/Maltab  Environment 

The  computation  of  the  linkage  fluxes  ( Aqs  j  and  ( Ads  j  is  achieved  by  means  of  the  integration  of  the  quantities 
( v9*  -  rsiqs )  and  (vds  -  rsids ) ,  respectively.  An  accurate  integration  of  a  wavefonn  requires  the  exact  determination  of  the  initial 

conditions.  These  initial  conditions  should  be  determined  in  the  block  of  integration.  Otherwise,  a  constant  value  is  added  to  the  output 
wavefonn  of  the  integration  and  the  computed  torque  wavefonn  has  big  ripples.  The  desired  initial  conditions  of  the  magnetic  flux  are 
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not  easily  computed.  Nevertheless,  the  two  cascaded,  first  order,  low  pass  filters  at  the 
output  of  the  estimator  eliminate  the  ripples  and  improve  the  integration  procedure. 

The  final  step  of  any  design  is  the  establishment  of  its  quality,  as  well  as  the 
performance  and  the  reliability  of  it.  The  experimental  test  is  a  procedure  able  to  check 
whether  the  above  criteria  are  satisfied.  Using  the  state-of-the-art  hardware  and  software 
developed  in  the  power  electronics  laboratory  at  the  Naval  Postgraduate  School,  realistic 
conclusions  about  the  performance  of  the  electromagnetic  torque  estimator  can  be  drawn. 
A  schematic  description  of  the  laboratory  experiment  is  demonstrated  in  Figure  2. 


Figure  2.  Schematic  Description  of  the  Laboratory  Experiment 


The  induction  motor  is  connected  through  the  proper  power  electronics  with  the 
Field  Programmable  Gate  Array.  The  electromagnetic  torque  estimator  model  is 
converted  into  XILINX  blocks  using  the  XILINX  library  contained  in  the 
Simulink/Matlab  library  browser.  By  means  of  the  appropriate  XILINX  software,  the 
Field  Programmable  Gate  Array  is  programmed  to  execute  the  transfonnation  of  the 
signals  into  the  stationary  reference  frame  as  well  as  the  electromagnetic  torque 
calculation.  Consequently,  as  the  induction  motor  accelerates  or  decelerates,  the  FPGA 
calculates  the  electromagnetic  torque  continuously.  The  results  are  plotted  and  analyzed 
on  the  computer. 
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The  experimental  application  poses  certain  problems  caused  by  the  external  noise 
interference  and  the  appearance  of  constant  offsets  at  the  waveforms.  The  main  problem, 
though,  is  related  to  the  numerical  integration.  The  flux  is  estimated  through  the 
integration  of  the  stationary  voltage  and  the  current  wavefonns.  However,  if  there  is  an 
offset  at  the  input  of  the  integrator,  then  a  ramp  error  occurs  at  the  output  of  the 

integration.  The  application  of  the  operator  instead  of  the  operator  (J/)  gives  a 

stable  integration  output.  The  flux  integrator  in  XILINX  software  is  thus  going  to  be 
implemented  as  ( Furthermore,  the  operator  1  is  transformed  into  the  digital 


1-z 

domain  using  the  Euler  approximation  ( s  = - )  and  it  is  implemented  with 

elementary  operations  such  as  the  time  delay,  the  sum  and  the  scaling. 

The  electromagnetic  torque  estimator  in  XILINX  environment  is  shown  in  Figure 
3: 


xxi 


integrator  1/(sM) 


Delay  1 


Figure  3.  Torque  Estimator  in  XILINX  Environment 

The  accuracy  of  the  electromagnetic  torque  estimator  is  verified  by  an  electrodynamometer.  The  electrodynamometer  is 
connected  to  the  motor  through  a  belt.  It  applies  a  mechanical  torque  to  the  motor,  which,  in  ideal  conditions,  is  equal  to  the 
electromagnetic  torque.  In  practice,  the  measured  electromagnetic  torque  should  be  a  little  bit  larger  in  order  to  surpass  the  mechanical 
load  and  the  losses.  Since  the  difference  is  insignificant,  though,  the  electrodynamometer  is  a  good  way  to  test  the  electromagnetic 
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torque  estimator.  A  zero  load  condition  has  an  initial  electromagnetic  torque  of 
(lAlbf.in  ).  This  value  must  certainly  be  added  to  the  applied  mechanical  load.  The  final 
value  deriving  from  this  addition  should  thus  be  compared  with  the  estimator 
measurement. 


Figure  4.  Electromagnetic  Torque  Response  after  Applying  Three  Different  Step  Loads 

Three  different  step  loads  are  applied  at  the  motor,  as  seen  in  Figure  4.  The 
calculated  steady  state  torque  approximates  the  electrodynamometer  settings.  The 
analysis  of  the  results  shows  that  the  torque  estimator  works  successfully  for  the 
induction  motor. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Induction  machines  are  used  worldwide  in  many  residential,  commercial, 
industrial  and  utility  applications.  They  transform  electrical  energy  into  mechanical 
energy.  An  induction  motor  may  be  part  of  a  pump  or  fan,  or  connected  to  some  other 
form  of  mechanical  equipment  such  as  a  winder,  a  conveyor  belt  or  a  mixer.  Induction 
motors  have  existed  for  many  years  but  were  always  limited  in  their  application  because 
it  was  difficult  to  control  the  speed  of  the  motor.  Nowadays  induction  motors  are  the 
preferred  choice  among  the  industrial  motors  due  to  their  rugged  construction,  the 
absence  of  brushes  (which  are  required  in  most  DC  motors),  the  modern  power 
electronics,  and  the  ability  to  control  the  speed  of  the  motor. 

B.  OBJECTIVE 

There  are  certain  kinds  of  applications  that  require  the  design  of  a  very  efficient 
estimator  of  the  induction  motor  electromagnetic  torque.  An  application  of  this  kind  is  the 
electric  vehicle.  Knowing  the  electromagnetic  torque  we  are  able  to  control  it. 
Furthermore,  we  can  control  the  speed  of  the  vehicle  faster  and  more  stably  after  adding  a 
torque  regulator  to  the  speed  loop.  The  accelerator  of  the  electric  vehicle  could  be  a 
closed  loop  speed  control  system  that  would  compute  the  electromagnetic  torque,  and 
would  accordingly  use  that  computation  in  order  to  control  the  torque  and  thus  accelerate 
and  decelerate  the  vehicle. 

The  objective  of  this  study  is  to  design  an  accurate  electromagnetic  torque 
estimator  of  a  three-phase  induction  motor,  without  the  use  of  any  sensor.  The  designed 
torque  estimator  without  any  sensor  is  based  on  the  mathematical  model  of  the  motor.  It 
can  easily  be  used  to  any  control  scheme  without  any  influence  on  the  efficiency  of  the 
motor. 
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c. 


INDUCTION  MOTOR  OVERVIEW 


An  induction  motor  is  an  electric  motor  that  is  driven  by  an  alternating  current.  It 
consists  of  two  basic  parts:  an  outside  stator  having  windings  supplied  with  alternating 
current  to  produce  the  rotating  magnetic  field,  and  an  inside  rotor  attached  to  the  output 
shaft,  as  shown  in  Figure  5. 


Figure  5.  A  Three-phase,  wye-connected  Induction  Motor  Schematic  and  a  Real 

Induction  Motor 

There  are  two  types  of  induction  machine  rotors:  squirrel  cage  rotors  and  wound 
rotors.  In  this  application  a  three-phases,  four  poles  squirrel  cage  rotor  induction  motor  is 
going  to  be  used.  In  this  kind  of  rotor  the  current  flows  through  copper  or  aluminum  bars. 
The  bars,  which  are  uniformly  distributed  and  embedded  in  a  ferromagnetic  material,  all 
end  up  at  a  common  ring  at  each  end  of  the  rotor,  as  shown  in  Figure  6.  This  rotor 
structure  is  very  simple,  robust  and  reliable. 


Figure  6.  Squirrel  Cage  Rotor 
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As  a  consequence  of  the  common  ring  existence  at  each  end  of  the  rotor,  these 
aluminum  conductors  are  short  circuited.  The  current  can  thus  flow  up  and  down  these 
aluminum  bars.  The  operating  principle  of  the  induction  motor  is  based  on  these  currents. 
The  stator  windings  create  a  rotating  magnetic  field.  This  changing  magnetic  field  causes 
flux  through  these  short-circuited  loops  of  the  rotor  and  induces  current.  The  induced 
currents  create  another  magnetic  field  that  interacts  with  the  magnetic  field  created  by  the 
stator  and  in  effect  causes  a  rotational  motion  in  the  rotor  [1], 

As  a  result  of  this  operating  principle  the  rotor  cannot  rotate  at  the  same  speed  as 
the  stator  rotating  field.  If  that  happens,  there  will  be  neither  change  in  the  flux  nor 
induced  voltage,  and  consequently  no  induced  currents  in  the  rotor.  Therefore,  there  will 
not  be  any  reaction  force  against  the  stator  rotating  magnet  motive  force.  Consequently, 
the  rotor  must  always  rotate  more  slowly  than  the  stator  rotating  field,  namely  “slip.”  The 
question  of  how  fast  the  rotor  slips  depends  on  the  balance  between  the  electromagnetic 
torque  of  the  motor  and  the  torque  load  [2]. 

D.  APPROACH 

An  accurate  estimator  of  the  electromagnetic  torque  of  a  three-phase  symmetrical 
induction  motor  is  going  to  be  designed  without  the  use  of  any  sensor.  This  goal  is  going 
to  be  achieved  by  building  the  mathematical  model  of  the  induction  motor  and  deriving 
the  appropriate  mathematical  equations  that  describe  it.  Using  the  proper  transfonnations 
to  reduce  the  complexity  of  these  differential  equations,  the  electromagnetic  torque 
equation  of  the  induction  motor  is  going  to  be  derived.  The  designed  estimator,  using 
Simulink/Matlab,  will  be  based  on  that  equation.  The  computer  simulation  will  detennine 
the  transient  and  steady  state  torque  of  the  motor.  Furthermore,  the  computation  of  the 
torque  of  the  real  motor  and  the  proof  of  the  efficiency  of  our  estimator  will  be  done  by 
converting  the  model  in  XILINX  environment  and  using  the  XILINX  FPGA  (Field 
Programmable  Gate  Array). 

The  thesis  is  organized  as  follows.  After  the  introduction  of  Chapter  I,  Chapter  II 
contains  the  derivation  of  the  torque  equation,  computer  simulation  and  simulation 
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results.  Chapter  III  describes  the  real  implementation  of  the  torque  estimator  and  the 
hardware.  Moreover,  Chapter  III  analyzes  the  waveforms  that  have  derived  from  the  real 
application.  Finally,  Chapter  IV  includes  general  conclusions  and  possibilities  for  further 
research. 
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II.  OBSERVER  OF  THE  ELECTROMAGNETIC  TORQUE 


A.  DERIVATION  OF  TORQUE  EQUATION 

1.  Flux  Linkage  and  Voltage  Equations  in  the  Arbitrary  (qdO)  Reference 
Frame 


Stator  Windings  .  7  ^  n  ^ 

Equivalent  Rotor  W ladings 


Figure  7.  Circuitry  Representation  of  the  Squirrel  Cage  Induction  Motor 

A  circuitry  representation  of  the  squirrel  cage  induction  motor  is  shown  in  Figure 
7.  It  is  convenient  for  the  analysis  of  the  circuit  to  transform  all  rotor  variables  into  the 
stator  windings  by  appropriate  turn  ratios  [1]. 


The  voltage  equations  become: 
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In  our  model,  the  total  flux  linking  each  of  the  three  stator  phases,  as  well  as  the  rotor  phases,  is  a  combination  of  the  six 
currents.  Consequently,  the  flux  linking  each  stator  phase  winding  can  be  written  as  [1]: 
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In  the  above  equations,  Lh  and  Lms  respectively  stand  for  the  leakage  and  magnetizing  inductances  of  the  stator  windings; 
Llr  and  Lmr  are  defined  as  the  leakage  and  magnetizing  inductances  for  the  rotor  windings,  and  the  inductance  Lsr  is  the  amplitude  of 
the  mutual  inductances  between  the  stator  and  the  rotor  windings.  The  quantity  9r  is  the  electrical  angular  displacement  of  the  rotor 
with  respect  to  the  (as)  axis. 
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We  can  combine  the  flux  linkages  of  the  stator  in  a  matrix  form  as  follows: 

(10) 
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The  winding  inductances  in  matrix  form  are  [1]: 
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The  complexity  of  our  mathematical  analysis,  which  is  required  for  the  derivation 
of  the  electromagnetic  torque  equation,  can  be  reduced  by  transfonning  the  variables 
associated  with  the  stator  and  the  rotor  elements  into  the  arbitrary  (qdO)  reference  frame. 
In  general,  the  desired  variables  can  be  transfonned  into  the  arbitrary  (qdO)  reference 
frame  by  multiplying  them  by  the  proper  transfonnation  matrices. 
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The  proper  transformation  matrices  for  the  stationary  and  rotor  variables  are 
respectively  [1]: 
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In  the  above  transfonnation  matrices  9  is  the  angular  displacement  of  the  arbitrary 
reference  frame. 

A  change  of  variables  that  formulates  a  transformation  of  variables  of  stationary 
circuit  elements  to  the  arbitrary  reference  frame  can  be  expressed  as: 


(16) 
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where  f  is  an  arbitrary  variable. 

The  following  result  is  derived 
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This  can  be  seen  from  equation  (10) 
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Multiplying  both  sides  by  K  and  transforming  the  stator  and  rotor  currents  into 
the  arbitrary  reference  frame  using  equation  (16),  it  is  obtained 
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The  left  hand  side  is  the  flux  linkages  of  the  stator  in  the  arbitrary  (qdO)  reference 
frame  in  a  matrix  form.  Knowing  that  KsLsK~'  and  KsLsrK~'  are  diagonal  it  is  obtained 
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and  equation  (17)  follows  easily. 

The  detennination  of  the  stator  voltage  expressions  in  the  arbitrary  (qdO) 
reference  frame  is  required  for  the  derivation  of  the  electromagnetic  torque  equation. 

Equation  (2)  can  be  represented  in  a  matrix  form  as  follows 
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Multiplying  both  sides  by  Ks  and  transforming  the  stator  currents  and  the  linkage 
fluxes  into  the  arbitrary  (qdO)  reference  frame  it  is  obtained: 
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Computing  the  derivative  at  the  right  hand  side,  it  is  derived 
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The  left  hand  side  represents  the  stator  voltages  in  the  arbitrary  reference  frame  in 
a  matrix  form.  Knowing  that  KsKj  is  the  unity  matrix  and  computing  the  term 
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and  expanding  the  matrices  it  follows  that  the  stator  voltage  equations  in  the  arbitrary 
reference  frame  are 
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2.  Electromagnetic  Torque  Equation 

The  total  energy  stored  in  a  rotating  electromagnetic  system  with  J  electrical 
inputs  can  be  expressed  as  [1]: 
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»7(<’„<2.-.^)=iESi«v,  <19> 


z  p=i  ?=i 


In  the  case  of  the  three-phase  induction  motor,  there  exist  two  electrical  inputs; 
the  stator  and  the  rotor.  Equation  (19)  is  reduced  as  [1]: 

1  1 
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The  tenn  Zn  describes  the  partial  derivative  of  the  linkage  flux  of  the  first  input 
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the  motor  only  is  derived  by  combining  the  equations  (4),  (7)  and  (20): 
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Following  the  same  procedure  through  all  three  phases  of  the  induction  motor,  the 
total  energy  stored  in  the  coupling  fields  can  be  written  as: 

i  \L  i.  +-L  (7  2 +4  2 +4  2)  (21) 

cs  J  sr  br  ^  r  y  ar  br  cr  J  V  / 

_cj 

In  rotational  systems,  the  electromagnetic  torque  is  defined  as  the  derivative  of 
the  co-energy  IT,  with  respect  to  the  actual  angular  displacement  of  the  rotor 0rm. 

Assuming  that  the  motor  is  magnetically  linear  then  the  field  energy  is  equal  to  the  co¬ 
energy.  Furthermore,  the  actual  angular  displacement  of  the  rotor  for  a  P-pole  motor  is 
related  to  the  electric  angular  displacement  as  [1]: 
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e  =  —  e  (22) 
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So,  the  electromagnetic  torque  may  be  expressed  as  [1]: 
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using  equation  (20)  to  obtain 
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Equation  (23)  can  be  transformed  into  the  arbitrary  reference  frame  as  [1]: 
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Solving  equation  (17)  to  find  i  and  idr  the  following  is  derived: 
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Substituting  equations  (25)  and  (26)  for  equation  (24),  the  final  simple  equation  of  the 
electromagnetic  torque,  which  is  a  combination  of  the  transformed  stator  currents  and 
linkage  fluxes,  is  derived: 


B.  TORQUE  ESTIMATOR  IN  SIMULINK/MATLAB  ENVIROMENT. 


1.  Designing  the  Torque  Estimator 


p  'l 

Equation  T  =  -~(iqsAds  -  idsAqs ) 


shows  that  the  electromagnetic  torque  depends 


on  the  stator  currents  and  the  linkage  fluxes.  The  stator  currents  are  measurable 
quantities.  There  is  electronic  equipment  that  measures  the  current  with  great  accuracy. 
On  the  other  hand,  the  linkage  flux  is  a  quantity  that  is  very  difficult  to  be  measured.  This 
observation  leads  to  the  conclusion  that  the  more  accurate  the  estimation  of  the  flux  the 
more  efficient  the  torque  calculation.  The  design  of  an  appropriate  sensor,  which  could 
measure  the  flux  and  could  be  placed  between  the  stator  and  the  rotor,  would  solve  the 
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problem.  However,  this  idea,  besides  the  constructional  problems,  could  affect  the 
magnetic  fluxes  in  the  motor  and  reduce  its  efficiency.  On  the  other  hand,  torque 
estimation  without  any  sensor,  based  on  the  mathematical  model  of  the  motor,  could 
easily  be  used  without  any  influence  on  the  flux. 


The  analysis  of  the  motor  and  the  derivation  of  the  voltage  and  magnetic  flux 
equation  have  been  done  to  the  arbitrary  reference  frame.  Nevertheless,  a  closer 
examination  of  voltage  equations  (18)  shows  that  their  transformation  from  the  arbitrary 
into  the  stationary  reference  frame  renders  the  desired  linkage  flux  computation  easier. 
The  speed  of  the  stationary  reference  frame  is  zero.  In  other  words,  the  voltage  equations 
are  more  easily  solved  by  replacing  (0)  =  £ ,JL)  with  zero.  The  q-axis  and  d-axis  linkage 

dt 

fluxes  in  the  stationary  reference  frame  are  respectively  [5]: 
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which  are  derived  by  solving  equations 
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The  torque  equation  in  the  stationary  reference  frame  becomes: 

T  =  ~(iqsAds  -idsAqs)  =  ~{iqs  (J (v*  -rsids)dt^-ids  (J  (vqs  -rf  s)dt^  (30) 
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The  final  torque  equation  depends  on  the  stator  voltage,  stator  currents  and  the 
resistance  of  stator  phases — quantities  that  can  be  measured  accurately.  Actually,  the 
variation  of  the  stator  resistance  with  time/temperature  is  significant.  The  equations  are 
very  sensitive  to  the  resistance,  especially  at  low  speed. 

Accordingly,  the  torque  estimator  is  fonnulated  based  on  equation  (30)  and  in  the 
Simulink/Matlab  environment  it  can  be  designed  as  in  Figure  8: 


Electromagnetic  Torque 


Figure  8.  Torque  Estimator  in  Simulink/Matlab  Environment 

Standard  Simulink  blocks  are  used  for  the  implementation  of  this  model.  It 
consists  of  four  inputs,  the  stator  voltage  and  the  stator  currents  in  the  stationary 
reference  frame,  and  the  desired  electromagnetic  torque  as  its  output. 

2.  Computer  Simulation — Analysis  of  the  Results 

A  computer  simulation  model  is  going  to  be  used  in  order  to  test  the  estimator. 
This  model  simulates  a  four  pole,  60  Hz,  three-phase  induction  motor  with  the  following 
rating  parameters: 
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Power 

Phase  voltage 

tpm 

rs 

rr 

4 

L 

ms 

4 

3  hp 

220  Volts 

1710 

0.435 

ohms 

0.816 

ohms 

0.002  H 

0.0462  H 

0.002  H 

Table  1.  Parameters  of  the  Computer  Simulated  Three-phase  Induction 

Motor 


The  computer  simulation  model  was  designed  in  a  way  that  allows  us  to  choose 
the  appropriate  frequency  speed  and  consequently  to  define  the  reference  frame. 
Therefore,  the  four  desired  inputs  to  the  estimator  ( vqs ,  vds ,  iqs ,  ids ) ,  which  are  extracted 

from  the  computer  model,  are  already  transformed  into  the  stationary  reference  frame  [4]. 

The  induction  motor  is  simulated  for  a  time  (tstop  =  2  sec)  with  the  sample  time 
of  (tsample  =  0.00002  sec).  A  step  load  torque  is  applied  at  (tstep  =  0.5  sec)  on  the  motor. 
Ideally,  in  the  case  of  a  motor  driving  a  load,  the  load  torque  is  equal  to  the 
electromagnetic  torque.  Consequently,  the  applied  load  torque  affects  the  electromagnetic 
torque.  The  results  are  displayed  and  analyzed  on  a  scope. 


<U 


3 


s 

SP 


s 


Figure  9.  Electromagnetic  Torque  after  Applying  a  Step  Load  Torque  on  the  Motor 
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The  output  of  the  designed  torque  estimator,  as  shown  in  Figure  9,  is  not  the 
expected  output.  After  the  application  of  the  step  function,  in  steady  state,  it  should  be 
constant  without  the  ripples  that  appear  in  our  case.  Furthermore,  in  computer  simulation 
the  waveforms  are  perfect  and,  most  importantly,  they  are  not  disturbed  by  the 
interference  of  any  undesired  external  noise  that  could  create  these  ripples. 

The  undesired  ripples  at  the  output  are  caused  by  the  integration.  The  computation 
of  the  linkage  fluxes  ( Aqs )  and  (Ad' )  is  achieved  by  the  integration  of  the  quantities 

(vqs  -  r/s )  and  (v*  -  rf )  respectively.  An  accurate  integration  of  a  waveform  requires 
the  exact  determination  of  the  initial  conditions. 

Assuming  that  a  sine  waveform  is  integrated,  it  can  be  written  [6]: 


Initial  Conditions  j 


These  initial  conditions  should  be  determined  in  the  block  of  integration. 
Otherwise,  a  constant  value  is  added  to  the  output  waveform  of  the  integration. 


X:  0.1207 
Y:  0.4683 


q-axis  linkage  flnx  ls  present 


A  small  offset  of  -0.0073 


-0.5 


0.1 


0.11  0.12  0.13  0.14 


0.15  0.16 


An  offset  of  -0.476  i: 


-1 


0.1  0.11  0.12  0.13  0.14  0.15  0.16 

Time  in  seconds 


Figure  10.  The  Waveforms  of  the  Linkage  Fluxes  Shifted  Vertically  by  Constant  Offsets 
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The  subtraction  of  the  two  waveforms  with  offset  causes  the  ripples  that  appear  in 
the  output  torque  of  Figure  10.  The  desired  initial  conditions  of  the  magnetic  flux  are  not 
easily  computed.  Instead,  there  are  two  possible  solutions  against  the  integration  problem. 

a.  Applying  a  High  Pass  Filter  after  the  Integration 

The  application  of  a  first  order  high  pass  filter  before  the  integration 
blocks  eliminates  the  offset  values  from  the  flux’s  waveforms  and  hence  improves  the 
integration  procedure.  The  fundamental  frequency  of  the  waveforms  is  60  Hz.  So  a  high 
pass  filter  with  a  cutoff  frequency  of  10  Hz  is  selected  as  shown  in  Figure  11. 


Figure  1 1 .  Applying  a  high  pass  filter  after  the  integrator 

The  magnitude  of  the  frequency  response  of  a  first  order  high  pass  filter 
$ 

with  transfer  function  (Tf  =  — - — )  can  be  displayed  via  a  bode  plot  as  shown  in  Figure 

5  +  10 


12. 
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Figure  12.  Frequency  Response  of  a  Fligh  Pass  Filter  with  Transfer  Function 

V  =  ’/,  + 10 

The  offset  of  the  flux  linkage  waveforms  has  been  removed  by  the  high 
pass  filters  as  shown  in  Figure  13. 


0.7  0.71  0.72  0.73  0.74  0.75  0.76 


Figure  13.  The  Effect  of  the  High  Pass  Filters  on  the  Flux  Linkage  Waveforms  in  Steady 

State 
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After  the  adjustment,  which  has  been  made  by  the  high  pass  filters,  the 
electromagnetic  torque  is  shown  in  Figure  14. 


Figure  14.  Electromagnetic  Torque  of  the  Motor  after  the  Placement  of  Two  High  Pass 
Filters  before  the  Integration  Blocks  on  the  Torque  Estimator  Model 

b.  Placing  a  Low  Pass  Filter  in  the  Output  of  the  Torque  Estimator 

The  analysis  of  the  undesired  waveform  of  Figure  9  in  the  frequency 
domain  can  be  made  by  computing  the  discrete  Fourier  transform  of  the  signal.  The 
discrete  Fourier  transform  is  an  algorithm  used  in  situations  where  there  is  not  an 
analytical  expression  of  a  signal  and  hence  its  frequency  spectrum  needs  to  be 
determined.  An  efficient  way  of  computing  the  discrete  Fourier  transform  is  the  fast 
Fourier  transform.  By  applying  the  fast  Fourier  transform  on  the  signal  of  Figure  9,  the 
spectrum  in  Figure  15  is  obtained. 
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Figure  15.  Fast  Fourier  Transform  of  the  Signal  Shown  in  Figure  9 

The  analysis  of  the  frequency  spectrum  illustrated  in  Figure  15  shows  that 
the  signal  contains  a  constant  value  quantity  at  0  Hz  and  has  a  fundamental  frequency  of 
60  Hz.  In  order  to  extract  the  constant  quantity  that  in  fact  represents  the  electromagnetic 
torque,  a  low  pass  fdter  can  be  applied  to  the  output  of  the  torque  estimator.  As  the 
fundamental  frequency  is  at  60  Hz,  a  second  order  low  pass  filter  with  cutoff  frequency 
of  10  Hz  can  be  used  to  eliminate  the  unvarying  quantity.  The  transfer  function  of  the  low 
pass  filter  can  be: 

Tf-  ,  100 

52 +20^  +  100 

The  bode  plot  of  the  second  order  low  pass  filter  is  shown  in  Figure  16. 
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Figure  16.  Frequency  Response  of  a  Second  Order  Filter  with  Transfer  Function 

7V  =  100/ 

J  A2 +205 +  100 


The  electromagnetic  torque  after  the  low  pass  filter  becomes  as  shown  in 

Figure  17. 


Figure  17.  Electromagnetic  Torque  after  Applying  a  Second  Order  Low  Pass  Filter  to  the 

Output  of  the  Estimator 
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c. 


Comparing  the  Two  Proposed  Solutions 


Figure  18.  Comparing  the  Torque  Diagrams  after  Applying  a  Low  Pass  Filter  to  the 
Output  of  the  Estimator  or  Two  High  Pass  Filters  before  the  Integration  Blocks 

The  implementation  of  two  high  pass  filters  affects  the  response  of  the 
electromagnetic  torque  as  it  reaches  the  steady  state  value  faster,  as  seen  in  Figure  18.  On 
the  other  hand,  the  steady  state  values  are  almost  the  same,  proving  that  both  solutions 
can  be  used  with  accurate  steady  state  results.  The  implementation  of  the  estimator  in  the 
computation  of  the  electromagnetic  torque  of  a  real  induction  motor  requires  that  it 
should  be  as  simple  as  possible.  It  is  easier  to  construct  a  low  pass  filter  compared  to  a 
high  pass  filer.  As  a  result,  the  final  modification  of  the  torque  calculator  is  composed  of 
an  extra  second  order  low  pass  filter  at  the  end  of  it.  The  accuracy  of  the  designed 
estimator  is  analyzed  in  the  next  chapter.  Using  the  appropriate  hardware  and  software, 
the  torque  of  a  real  three-phase  induction  motor  is  computed  continuously  as  the  motor 
runs,  giving  an  estimation  of  the  precision  of  our  torque  calculator. 
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III.  COMPUTING  THE  TORQUE  OF  A  REAL  MOTOR 


A.  HARDWARE  DESCRIPTION 

The  final  step  of  any  design  is  the  establishment  of  its  quality,  and  the 
performance  and  the  reliability  of  it.  The  experimental  test  is  a  procedure  able  to  check 
whether  the  above  criteria  are  satisfied.  Using  the  state  of  the  art  hardware  and  software 
developed  in  the  power  electronics  laboratory  at  the  Naval  Postgraduate  School,  realistic 
conclusions  about  the  perfonnance  of  the  electromagnetic  torque  estimator  can  be  drawn. 
A  description  of  the  laboratory  experiment  tools  is  demonstrated  below. 

1.  Four  Pole  Squirrel  Cage  Induction  Motor 

An  example  of  a  four  pole  squirrel  cage  induction  motor  is  shown  below  in  Figure 
19. 
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Figure  19.  Model  8821  Four  Pole  Squirrel  Cage  Induction  Motor 


The  characteristics  of  the  experimental  induction  motor  are  displayed  in  Table  2. 


Output 

Stator 

Frequency 

Full  load 

Full  load 

Dimensions 

Net 

Power 

Voltage 

speed 

current 

weight 

175  W 

120/208  V 

60  Hz 

1 670r/min 

1.2  A 

308  x  291  x440  mm 

13.5  Kg 

Table  2.  Characteristics  of  the  Laboratory  Induction  Motor  (Rating  Values) 
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2.  Electrodynamometer 


In  Figure  20,  the  electrodynamometer  is  connected  to  the  motor  via  a  belt.  It 
applies  mechanical  torque  to  the  motor  that  ideally  is  equal  to  the  electromagnetic  torque. 
Comparing  the  output  of  the  estimator  with  the  applied  torque  by  the  dynamometer  is  a 
way  to  verify  its  accuracy. 


Figure  20.  Electrodynamometer  Connected  to  the  Motor  Via  a  Belt 

3  Field  Programmable  Gate  Array  (FPGA) 

The  XILINX  Virtex-4TM  development  board  shown  in  Figure  21  incorporates  the 
(XS4VLX25-10SF363)  Field  Programmable  Gate  Array  (FPGA)  as  well  as  all  the 
essential  microelectronics  in  order  to  cooperate  with  multi-input  logic  blocks,  flip-flops, 
multiplexors,  memory  and,  finally,  the  clock  generation  section. 
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Figure  2 1 .  XILINX  Virtex-4™  development  board 


The  Field  Programmable  Gate  Array  (FPGA)  is  the  fundamental  electronic  tool 
used  in  our  experiment.  It  is  a  semiconductor  device  containing  programmable  logic 
components  and  programmable  interconnects.  The  distinguishing  features  of  a  FPGA  are 
its  property  of  being  programmable  by  the  designer  and  its  ability  to  be  easily 
reprogrammed.  Using  the  appropriate  software,  this  electronic  device  can  be  programmed 
to  execute  complex  mathematical  functions  as  well  as  combinations  of  them.  The 
designed  electromagnetic  torque  estimator  is  a  combination  of  mathematical  operations. 
Therefore,  it  can  be  implemented  via  the  FPGA.  Furthermore,  the  FPGA  is  able  to 
execute  codes  in  a  parallel  and  consequently  very  fast  manner.  It  can  thus  be  connected  in 
parallel  with  the  motor  as  it  runs  and  operates  simultaneously. 

4.  System  Interface  Box 

The  system  interface  box  or  the  FPGA  development  kit  box  in  Figure  22  contains 
the  following  electronic  tools: 
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•  A  power  supply  source  that  converts  the  input  power  from  AC  (alternate 
current)  to  DC  (direct  current)  that  is  important  for  the  safe  and  exact 
operation  of  the  boards. 

•  The  analogue  signal  interface  board,  which  consists  of  eight  input 
channels:  four  for  the  voltage  and  four  for  the  current.  Its  main  goal  is  to 
reduce  the  values  of  the  input  analogue  signals  by  the  appropriate  ratio 
before  their  conversion  from  analogue  to  digital. 

•  The  FPGA  interface  board,  which  contains  two  analogues  to  digital 
converters  that  sample  the  input  signals  for  digital  signal  processing 
operations.  Many  channels  are  also  included  which  can  be  used  as  outputs 
on  the  FPGA. 

•  The  JTAG  cable  is  used  to  program  and  reprogram  the  FPGA.  In  more 
general  terms,  the  communication  between  the  field  gate  array  and  the 
computer  occurs  through  this  specific  cable. 


FPGA  interface 
board 

— 


Xilinx  Virtex-4™ 
development  board 


Figure  22.  System  Interface  Box 
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5. 


Final  Experiment  Setup 


r. 


at 


^  ^  - 


System  interface 
box 


Motor  connected 
with  the 

electrodynamometer 


Figure  23.  The  Experiment  Setup  Station 


Induction  motor 


Figure  24.  Schematic  Diagram  of  the  Experiment  Arrangement 


30 


The  experiment  station  is  shown  in  Figure23  and  the  concept  of  the  experiment 
arrangement  is  illustrated  in  the  schematic  diagram  in  Figure  24.  Four  channels  of  the 
analogue  signal  interface  board  are  used,  two  for  the  first  and  second  phase  currents 
(4  and  4),  and  two  for  the  line-to-line  voltages  ( vab  and  vhc).  After  their  conversion  to 

digital  quantities,  these  four  signals  enter  the  FPGA.  The  Field  Programmable  Gate  Array 
has  been  programmed,  by  using  the  appropriate  computer  software,  to  execute  the 
transformation  of  the  signals  into  the  stationary  reference  frame  as  well  as  the 
electromagnetic  torque  calculation.  Consequently,  as  the  induction  motor  accelerates  or 
decelerates,  the  FPGA  calculates  the  electromagnetic  torque  continuously. 

XILINX,  a  leading  manufacturer  of  FPGAs,  has  developed  the  appropriate 
software  needed  for  the  programming  of  them.  The  proper  software  and  the  steps  that  are 
required  for  the  final  programming  of  the  transistor  are  depicted  in  the  Figure  25. 


Figure  25.  Programming  the  FPGA  Using  Computer  Based  Software 

The  electromagnetic  torque  estimator  model  needs  to  be  converted  in  XILINX 
blocks  using  the  XILINX  library  contained  in  the  Simulink/Matlab  library  browser.  Once 
the  model  has  converted,  the  System  Generator  software  produces  the  VHDL  code  (Very 
High  Speed  Integrated  Circuit  Hardware  Description  Language).  After  the  VHDL  is 
generated  suitable  for  the  FPGA  files  are  composed  by  the  ISE  foundation  software. 
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Through  the  JTAG  cable  the  final  project  is  uploaded  in  the  FPGA.  The  results  are 
displayed  on  the  computer  screen  through  the  ChipScope  Pro  analyzer  software  where 
they  can  be  analyzed  and  plotted  [7]. 

B.  XILINX  SOFTWARE  ANALYSIS 

The  Simulink/Matlab  model  in  Figure  8  needs  to  be  designed  in  XILINX 
environment.  The  designing  of  the  XILINX  library  is  structured  in  a  similar  way  to  the 
Simulink/Matlab  library,  though  there  are  some  differences  in  the  block  parameters.  In 
fact,  the  XILINX  library  seems  to  be  less  extended  compared  to  the  Simulink/Matlab 
library.  In  addition,  its  capability  is  limited. 

1.  Integration  in  XILINX  Environment 

The  experimental  application  presents  certain  problems  caused  by  the  external 
noise  interference  and  the  appearance  of  constant  offsets  at  the  waveforms.  The  flux 
estimation  is  achieved  by  the  integration  of  the  stationary  voltage  and  the  current 
waveforms.  However,  if  there  is  an  offset  at  the  input  of  the  integrator,  a  ramp  error 
occurs  at  the  output  of  the  integration.  For  instance,  the  integration  of  a  sine  wave  with  an 
offset  results  in  [6]: 

J (sin(ruf) + DC )dt  =  J  sin(cot)dt  +  J DCdt  =  — -cos(ruf)  +  (DC )  t  +  Initial  Conditions  (32) 
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Sine  wave  Y=5+200sin(2*pi*20*t) 


Integration  output 


Figure  26.  Ramp  Error  at  the  Output  of  Integration  Due  to  an  Offset  in  the  Signal 

The  ramp  error  (DC)t  causes  unstable  integration  as  shown  in  Figure  26.  This 

integration  in  turn  causes  saturation  of  the  software  and  undesirable  effects  in  the 
experimental  application.  Accordingly,  besides  the  initial  condition  problem  that  arises 
with  the  low  pass  filter  at  the  output,  the  unstable  integration  problem  needs  to  be 
overcome. 

The  integration  output  can  become  stable  by  applying  a  unity  feedback  as  shown 
in  Figure  27. 
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The  transfer  function  of  the  unity  feedback  closed  loop  system  becomes: 


Tf  = 


C(s)  G(s) 


R(s )  1  +  G(s)H(s)  1  +  j/  s  +  l 


(33) 


As  derived  from  the  analysis  above,  the  application  of  the  operator  ( |  ) 
instead  of  the  operator  (/)  gives  a  stable  integration  output  as  shown  in  Figure  28. 


Thus,  the  flux  integrator  in  XILINX  software  is  going  to  be  implemented  as  ( j/  +  ^ ). 


Sine  wave  Y=5+200sin(2*pi*20*t) 
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Integiation  output  after  applying  the  operator  l/(s+l) 


III 
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Figure  28.  The  Operator  ( ^ )  Makes  the  Integration  Stable 


The  implementation  of  the  integral  operator  Ys^\  'n  rca'  '•'mc  application  is  yet 

another  problem  to  be  sorted  out.  In  general,  the  digital  implementation  includes  several 
hardware  limitations,  such  as  limited  memory,  finite  precision,  and  limited  speed 
execution.  Operations  that  require  a  finite  amount  of  data  and  make  the  algorithm 
computable  are  necessary. 
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According  to  the  Euler  approximation  technique,  a  transfer  function  in  the 
differential  operator  (s)  can  be  transformed  into  a  discrete  time  transfer  function  in  the 
time  delay  operator  (z)  by  substituting  [4]: 


1  —  z  ' 


(34) 


In  equation  (34)  Ts  is  the  sampling  interval.  The  clock  frequency  of  the  FPGA  is  25  MHz 
and  consequently  the  sampling  interval  is  40x10  9  sec. 


By  substituting  equation  (34)  in  the  integral  operator 


it  is  obtained: 


1  _  1  1  _  4xl0~8  _  3.99xl08 

j+T-l-z_I  ~  l-z-'  1.00000004 -z1  “  l-0.99999996z  1 

- +  1  - 5-  +  1 

T  40x  10~9 


(35) 


In  time  domain  equation  (35)  is  implemented  by  the  difference  equation: 

y[n\  =  0.99999996y[n  -1]  + 3.99x10  8 x[n]  (36) 


Equation  (36)  yields  the  following  block  diagram  realization: 


Figure  29.  Block  diagram  realization  of  equation  (36) 
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The  block  diagram  above  in  Figure  29  represents  the  operator  +  \  with 

elementary  operations  such  as  the  time  delay,  sum  and  scaling,  which  use  the  minimum 
memory  storage.  It  can  be  used  in  our  real  time  application.  The  integrator  in  XILINX 
software  has  been  designed  based  on  this  concept  and  it  is  shown  in  Figure  30. 


integrator  1/(s+1) 


Figure  30.  Integrator  in  XILINX  Environment 

2.  Low  Pass  Filter  in  XILINX  Environment 

The  construction  of  the  low  pass  filter  at  the  output  of  the  electromagnetic  torque 
estimator  is  based  on  the  same  concept.  Instead  of  a  second  order  filter,  which 
complicates  the  design,  two  cascaded  first  order  low  pass  filters  with  transfer  function 

Tf  =  ^  are  implemented.  The  designed  filter  is  shown  in  Figure  3 1 . 

5  +  10 
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Delay  1 


Figure  3 1 .  Two  Cascaded  First  Order  Low  Pass  Filters  using  XILINX  Blocks 
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3. 


Final  Design  of  the  Torque  Estimator  in  XILINX  Environment 


Delay  1 


Figure  32.  Electromagnetic  Torque  Estimator  Converted  in  XIILINX  Environment 


The  complete  design  of  the  electromagnetic  torque  estimator  is  illustrated  in  Figure  32. 
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4.  Implementation  of  Transformations  using  XILINX  Software 


The  first  and  second  phase  currents  (4  and  4),  as  well  as  the  line-to-line  voltages 
(vafcandvic),  have  been  defined  as  the  four  inputs  on  the  analogue  interface  board. 

Nevertheless,  these  quantities  must  be  transformed  into  the  stationary  reference  frame. 
This  transformation  is  going  to  be  achieved  through  the  FPGA  by  designing  the  proper 
software  using  XILINX  blocks. 

The  equations  that  describe  the  q-axis  and  d-axis  stator  currents  as  a  function  of 
the  stator  currents  (ia  and 4)  are  respectively  [2]: 

2^3  .  71..  2V3  . 

4s  =  —  sm(0  +  -)ia  +  —  sin(0)i4  (37) 

2V3  tv..  2^3  ...  . 

4s  = — —  cos(6*  +  y)4  — —cos  (0)ib  (38) 

Substituting  (0  =  0)  in  equations  (37)  and  (38),  the  stator  currents  are  transfonned 
from  the  arbitrary  into  the  stationary  reference  frame.  The  XILINX  software,  which 
transforms  the  stator  currents  into  the  stationary  reference  frame,  is  based  on  equations 
(37)  and  (38)  is  shown  in  Figure  33. 


Negate 


Figure  33.  Software  that  Transforms  the  Stator  Currents  ( 4  and 4  )  into  the  Stationary 

Reference  Frame  (From  [2]) 
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The  equations  that  describe  the  q-axis  and  d-axis  stator  voltages  as  a  function  of 
the  stator  line-to-line  voltages  are  respectively  [2] : 

vqs  =  |  cos(0)vab  + |  sin(<9  +  ^ )vbc  (39) 

vds  =  \  sin(#K*  -  \  cos(^  +  (4°) 

3  3  6 

Consequently,  the  final  model,  which  transforms  the  line-to-line  stator  voltages 
into  the  stationary  reference  frame,  is  based  on  the  equations  (39)  and  (40)  as  seen  in 
Figure  34. 


SineCosine  1 


Figure  34.  Software  that  Transforms  the  Line-to-Line  Voltages  ( vah  and  vhc )  into  the 

Stationary  Reference  Frame  (From  [2]) 
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IV.  EXPERIMENTAL  RESULTS  AND  ANALYSIS 


The  analysis  of  the  experimental  results  can  be  done  on  the  computer  using 
ChipScope  Pro  software.  ChipScope  Pro  inserts  a  logic  analyzer  and  a  bus  analyzer  core 
directly  into  the  design.  This  allows  the  user  to  view  any  internal  signal  through  the 
programming  interface  and  analyze  it  with  the  ChipScope  Pro  Logic  Analyzer. 
Furthermore,  the  data  can  be  exported  to  Matlab  environment  for  measuring  and  plotting 

m. 

The  accuracy  of  the  electromagnetic  torque  estimator  is  going  to  be  verified  by 
the  electrodynamometer.  The  electrodynamometer  is  connected  to  the  motor  via  a  belt.  It 
applies  mechanical  torque  to  the  motor  that,  ideally,  is  equal  to  the  electromagnetic 
torque.  In  practice,  the  measured  electromagnetic  torque  should  be  a  little  bit  larger  in 
order  to  get  over  the  mechanical  load  and  the  losses.  However,  the  difference  is 
insignificant.  Thus,  the  electrodynamometer  is  a  good  way  to  test  the  electromagnetic 
torque  estimator.  On  the  electrodynamometer,  there  is  an  indicator  that  shows  the  applied 
mechanical  torque  in  pound  force  inches  ( lbf.in  ).  In  addition,  a  rotational  knob  allows 
the  user  to  apply  step  loads  to  the  motor. 

The  applied  stator  phase  voltage  is  V  hase  RMS  =  WO  Volt .  The  rating  value  is 

V phase  rms  =  120  Volt .  As  the  motor  is  connected  to  the  electrodynamometer,  even  if  there 

is  no  load  applied,  it  has  an  initial  electromagnetic  torque  in  order  to  turn  itself  as  well  as 
the  electrodynamometer.  This  torque  needs  to  be  measured.  Furthermore,  it  is  necessary 
for  this  value  to  be  added  to  the  applied  mechanical  load  and  this  addition  in  turn  should 
be  compared  with  the  estimator  measurement. 
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The  measured  initial  electromagnetic  torque  is  2Albf.in ,  as  seen  in  Figure35. 


Figure  35.  Measured  electromagnetic  torque  of  the  motor  in  zero  load  condition 

The  rating  stator  current  is  Istator  =1.2  Amperes .  According  to  this  limitation,  the 
maximum  mechanical  torque  that  can  be  applied  is  10  Ibf.in .  So  three  step  loads  below 
10  lbf.in  are  applied  to  the  motor  in  order  to  test  the  electromagnetic  torque  estimator. 

The  electromagnetic  torque  wavefonn  after  applying  a  3  Ibf.in  step  load  for  the 
time  interval  of  four  seconds  and  back  again  to  zero  load  is  displayed  in  Figure  36. 
Following  the  same  procedure,  the  torque  waveforms  after  applying  6/bf.in  and  10  lbf.in 
step  loads  are  shown  in  Figure  38  and  Figure  40  respectively. 

The  measured  steady  state  electromagnetic  torque  at  3  lbf.in ,  6  lbf.in  and 
10  Ibf.in  loads  is  displayed  in  Figure  37,  Figure  39  and  Figure  41  respectively.  The 
measured  values  approximate  the  corresponding  actual  values  and  that  verifies  the 
accuracy  of  the  designed  electromagnetic  torque  estimator. 
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Figure  36.  Electromagnetic  Torque  Waveform  after  Applying  3  Ibf.in  Load  and  Back 

again  to  Zero  Load 


Figure  37.  Transient  and  Steady  State  Response  of  the  Electromagnetic  Torque  after 

Applying  a  Step  Load  of  3  lbf.in 

As  derived  from  Figure  37  the  measured  steady  state  torque  is  5. 6 lbf.in, 
approximate  to  the  actual  value  of  5  A  lb f  .in  . 
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Figure  38.  Electromagnetic  Torque  Waveform  after  Applying  6  Ibf.in  Load  and  Back 

again  to  Zero  Load 


Ligure  39.  Transient  and  Steady  State  Response  of  the  Electromagnetic  Torque  after 

Applying  a  Step  Load  of  6  lbf.in 

The  measured  steady  state  torque  as  shown  in  Ligure  39  is  8.5  Ibf  .in  as  expected. 
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Figure  40.  Electromagnetic  Torque  Waveform  after  Applying  10  lbf.in  Load  and  Back 

again  to  Zero  Load 


Figure  41 .  Transient  and  Steady  State  Response  of  the  Electromagnetic  Torque  after 

Applying  a  Step  Load  of  \0  lbf.in 

As  derived  from  Figure  41  the  measured  steady  state  torque  is  12.5 lbf.in, 


approximate  to  the  actual  value  of  \2Albf  .in  . 
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Figure  42.  Response  of  the  Electromagnetic  Torque  to  Three  Different  Loads 

The  steady  state  values  of  the  three  different  step  loads  are  shown  in  Figure  42. 
The  calculated  steady  state  torque  approximates  the  electrodynamometer  settings.  The 
analysis  of  the  results  shows  that  the  torque  estimator  works  successfully  for  the 
induction  motor. 

A  comparison  of  the  measured  electromagnetic  torque  waveforms  with  the 
computer  simulation  results  cannot  be  accomplished,  due  to  the  unknown  machine 
parameters  and  inertia  of  the  laboratory  induction  motor. 
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V.  CONCLUSIONS  AND  SUGGESTIONS 


A.  CONCLUSIONS 

The  quality  and  reliability  of  any  design  are  proven  by  how  well  it  works  in  a  real 
application.  The  transition  from  an  ideal  case,  such  as  computer  simulation,  to  a  real 
application  includes  factors  that  may  have  undesirable  effects.  The  experiment  of  an 
electromagnetic  torque  estimator  using  a  real  induction  motor  is  a  real  application  which 
involves  undesirable  external  noise,  constant  offsets  at  the  waveforms,  as  well  as 
variation  in  the  temperature  of  the  windings  and  thus  in  the  stator  resistance. 
Nevertheless,  the  analysis  of  the  results  shows  that  the  torque  estimator  works  very  well 
for  an  induction  motor  at  the  frequency  of  60  Hz.  The  calculated  steady  state  torque  is 
very  close  to  the  electrodynamometer  settings.  In  conclusion,  at  the  speed  of  60  Hz  the 
variation  of  the  stator  resistance  is  not  dominant  and  does  not  affect  the  torque 
calculation.  Finally,  the  numerical  integration  works  successfully. 

The  designed  torque  estimator  can  be  used  in  any  application  that  requires  the 
calculation  of  the  torque  of  an  induction  motor  that  works  at  a  high  frequency.  In  fact,  the 
calculation  of  the  electromagnetic  torque  is  an  important  step  before  the  speed  control  of 
the  motor.  The  speed  of  the  motor  can  be  controlled  by  controlling  how  much 
electromagnetic  torque  is  applied  on  the  induction  motor.  Such  a  direct  torque  control 
scheme  works  very  successfully  only  if  the  estimation  of  the  torque  is  accurate. 

B.  SUGGESTIONS  FOR  FURTHER  RESEARCH 

At  low  frequencies,  the  estimation  of  the  flux  becomes  difficult.  Two  important 
negative  effects  render  the  designed  torque  estimator  less  efficient:  the  low  pass  filter, 
which  should  have  low  cutoff  frequency,  as  well  as  the  small  transition  region,  which 
causes  complexity  problems.  Therefore,  instead  of  a  low  pass  filter,  a  detector  could  be 
designed  to  detect  the  initial  conditions  of  the  flux  and  remove  them  before  integration. 
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An  equally  important  problem  caused  by  the  low  frequency  is  the  stator  resistance 
variation.  The  voltage  drop  across  the  stator  resistance  is  comparable  to  the  applied 
voltage  due  to  the  high  stator  current.  Therefore,  any  variation  in  the  stator  resistance 
would  cause  an  error  in  the  flux  estimation.  In  this  case,  an  algorithm  that  could  compute 
the  variations  of  the  resistance  would  improve  the  performance  of  the  electromagnetic 
torque  estimator. 
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APPENDIX:  COMPUTER  SIMULATION  MODEL 


Induction  machine  model 


Induction  machine  model/AC  source 
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Induction  machine  model/induction  Machine 


Induction  machine  model/induction  machine/Saturation 
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Induction  machine  model/Induction  machine/psi  Os 
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Induction  machine  model/Induction  machine/psi  me 


Induction  machine  model/induction  machine/psi  qr 


Induction  machine  model/induction  machine/psi  qs 


3  Product 


Induction  machine  model/speed 


electrical  radians  per  second 


Induction  machine  model/Torque  Ec 


omega  _r 


wbby  2  H 


motor 
Constant  1 


Switch 


poles  by  2J 


Initial  Matlab  file 


motor  =  0;  %when  motor  =  1  then  10  HP  machine,  motor  =  0  then  3  hp 
machine 

saturation  =  1;  %OFF  when  set  to  one,  ON  when  set  to  zero 
tfault_on  =  30.75/60;  %page  181  to  describe  fog.  4.13-1 
tf ault_of f  =  36.75/60; 

%omega  =  2*pi*60;  %Modify  this  variable  to  change  reference  frame 

omega  =  0;  %Modify  this  variable  to  change  reference  frame 

omega_b  =  2*pi*60; 

omega  in  =  2*pi*60; 

t_start  =  .1; 

tstep= . 000002 ; 

twopiby3  =  2*pi/3; 

wbby2H  =  omega_b/2 /0 . 5 ; 

poles  =  4; 

polesby2J  =  poles/2/ . 089; 

fault  or  loadchange  =0;  %Fault  when  1,  load  change  when  not=l 
if  fault  or  loadchange  ==  1 

fault  =  1;  %  set  to  large  value  to  disable  (like  100) 

load_change  =  100; 

else 

fault  =  100;  %  set  to  large  value  to  disable  (like  100) 

load_change  =  1; 

end 

if  motor  ==1 

%Parameters  from  section  4.11,  10  HP  machine 

tfault_on  =  100*30.75/60;  %page  181  to  describe  fig.  4.13-1 

tf ault_of f  =  100*36.75/60; 

tstep= . 00002 ; 

tstop= . 9; 

rs= .0453; 

rr  =  .0222; 

Xls  =.0775; 

Xm  =2 . 042; 

Xlr  =  .0322; 
rsbyXls  =  rs/Xls; 
rrbyXlr  =  rr/Xlr; 

Xaq  =  1 / ( 1 /Xm+1 /Xls+1 /Xlr ) ; 

Xad  =  Xaq; 

XaqbyXls  =  Xaq/Xls; 

XaqbyXlr  =  Xaq/Xlr; 

XadbyXls  =  Xad/Xls; 

XadbyXlr  =  Xad/Xlr; 

XaqbyXm  =  Xaq/Xm; 

XadbyXm  =  Xad/Xm; 

V_phase  =  1; 
psi_mqic=0 ; 
psi_mdic=0 ; 
psi_qsic=0 ; 
psi_dsic=0 ; 
psi_qric=0 ; 
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psi_dric=0 ; 
omegar_ic  =  0; 
elseif  motor  ==  0 

%Parameters  from  section  4.10,  3  HP  machine 

t_start  =  0; 

tstop=l . 6; 

rs= .435; 

rr  =  .816; 

Xls  =.754; 

Xm  =26.13; 

Xlr  =  .754; 
rsbyXls  =  rs/Xls; 
rrbyXlr  =  rr/Xlr; 

Xaq  =  1/ ( 1/Xm+l/Xls+l/Xlr) ; 

Xad  =  Xaq; 

XaqbyXls  =  Xaq/Xls; 

XaqbyXlr  =  Xaq/Xlr; 

XadbyXls  =  Xad/Xls; 

XadbyXlr  =  Xad/Xlr; 

XaqbyXm  =  Xaq/Xm; 

XadbyXm  =  Xad/Xm; 

V_phase  =  220*sqrt (2) /sqrt (3) ;  %Peak  value  used  in  Simulink  for  source 
model 

%Steady  state  analysis 

omega  rotor  =  .  999*omega  in; 

slip  =  (omega  in-omega  rotor) /omega  in; 

Z1  =  rs  +  j*omega  in/omega  b*Xls; 

Z2  =  j*omega  in/omega  b*Xm; 

Z3  =  rr/slip  +  j*omega  in/omega  b*Xlr; 
las  =  V  phase/sqrt (2) / ( Z 1  +  Z2*Z3/ (Z2  +  Z3) ) ; 

Iar  =  -las* (Z2*Z3/ (Z2+Z3) ) /Z3; 

Te  =  3*poles/2 *Xm/omega_b*real ( j *con j ( las ) * Iar ) ; 

Iqse  =  sqrt (2) *real (las) ; 

Idse  =  -sqrt (2) *imag (las) ; 

Iqre  =  sqrt (2 ) *real ( Iar) ; 

Idre  =  -sqrt (2 ) *imag ( Iar) ; 
psi  mqe  =  Xm* ( Iqse+Iqre) ; 
psi  mde  =  Xm* ( Idse+Idre) ; 
psi  qse  =  Iqse*Xls  +  psi  mqe; 

psi  dse  =  Idse*Xls  +  psi  mde; 

psi  qre  =  Iqre*Xlr  +  psi  mqe; 

psi  dre  =  Idre*Xlr  +  psi  mde; 

Te  =  psi_dse*Iqse  -  psi_qse*Idse; 

psi  mqic=psi  mqe; 

psi  mdic=psi  mde; 

psi_qsic=psi_qse; 

psi_dsic=psi__dse; 

psi_qric=psi_qre; 

psi_dric=psi_dre; 

omegar  ic  =  omega  rotor; 

end 
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